Quantum-mechanical description of in-medium fragmentation 
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We present a quantum-mechanical description of quark-hadron fragmentation in a nuclear envi- 
ronment. It employs the path-integral formulation of quantum mechanics, which takes care of all 
phases and interferences, and which contains all relevant time scales, like production, coherence, for- 
mation, etc. The cross section includes the probability of pre-hadron (colorless dipole) production 
both inside and outside the medium. Moreover, it also includes inside-outside production, which is 
a typical quantum- mechanical interference effect (like twin-slit electron propagation). We observe 
a substantial suppression caused by the medium, even if the pre-hadron is produced outside the 
medium and no energy loss is involved. This important source of suppression is missed in the usual 
energy-loss scenario interpreting the effect of jet quenching observed in heavy ion collisions. This 
may be one of the reasons of a too large gluon density, reported by such analyzes. 

PACS numbers: 24.85.-|-p, 12.40.Gg, 25.40.Ve, 25.80.Ls 



I. INTRODUCTION 

Hadronization in a nuclear environment has always 
been a precious source of information about the space 
time pattern of hadronization. This process is character- 
ized by the production length of a pre-hadron (a colorless 
dipole), with the subsequent development of the hadronic 
wave function. A perturbative description of this process 
as radiation of a qq pair and creation of colorless dipole 
evolving to the final pion, as is illustrated in Fig. 1, was 
proposed in [1, 2]. Although the model was partly prob- 




FIG. 1: Reaction j*p nX. The incoming virtual proton 
is absorbed by a valence quark of a bound nucleon, leading 
to the production of a quark and gluon, — > qig. The 
gluon decays to a 52 53, and gi — 52 fuse to the final pion. 
The intermediate and final partons can experience multiple 
interactions in the nuclear medium. 

abilistic, it led to rather successful description [2] and 
even prediction [1] for semi-inclusive hadron production 
in deep- inelastic scattering (DIS) [3]. A pure classical 
string model [4-6] fitted to data also leads to a good 



agreement. No fully quantum mechanical description of 
in- medium hadronization has been available so far. 

The widely debated question whether the hadroniza- 
tion process ends by a leading pre-hadron production 
"within or without the medium?" [2], strictly speaking 
does not have a definite answer. In quantum mechan- 
ics a pre-hadron may be created both inside and outside 
the medium, and interference of the corresponding am- 
plitudes is important, as is demonstrated below. 

Here we are developing the model [1, 2] of perturbative 
hadronization and employ the Berger model [7] , which we 
improved in a recent paper [8] . In this model an energetic 
quark produced in a hard reaction, like e^e~ annihila- 
tion, DIS, or high-pr scattering, creates a leading pion, 
which carries a major fraction ^ 1 of the quark mo- 
mentum, via perturbative radiation of a gluon decaying 
into a qq pair. Then the q fuses with the original quark 
into the pion, as is illustrated in Fig. 1 for 7*-nucleus 
collision. In general, fragmentation is process dependent, 
due to the higher twist terms calculated in [7, 8], which 
we neglect here. Therefore our results are applicable to 
any hard reaction leading to the production of a quark 
jet. 

The main approximations in the Berger model [7] are: 
(i) the calculation is done to lowest order of (1 — Zh), 
which is considered as a small parameter; (ii) the pion is 
treated as a qq state with no transverse or longitudinal 
motion of the quarks. New calculations without these 
approximations were done in [8] , where the cross section 
of the process 7*p ^ ttX was calculated to all orders of 
(1 — Zh) and with a realistic model for the light-cone pion 
wave function 

The calculations performed in [7, 8] were done in mo- 
mentum representation, employing the Feynman diagram 
technique, and for the case of a proton target. Usually 
for nuclear targets the impact parameter representation 
is more effective, since at high energies the impact pa- 
rameters do not vary during propagation through the 
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nucleus, which aUows to apply a Glauber- like eikonal- 
ization. If the energy is not sufficiently high to freeze the 
impact parameters, one should integrate over all possible 
paths of the propagating partons. 

Here we employ the light cone Greens function for- 
malism which is the essential tool for the calculation of 
the nuclear effects [9-12]. We calculate the ratio oi pt- 
integrated cross sections, 



RA/p{Zh) 



dai"f*A^ TTX)/dzh 
Ada{-f*N ^nX)/dzh' 



(1) 



as function of pion fractional light-cone momentum Zh- 

We consider gluon decay and pre-hadron (5291) produc- 
tion both inside and outside the nucleus. However, these 
two possibilities can be clearly separated only in a prob- 
abilistic approach. In fact, in quantum mechanics the 
cross section is related to the square of the process am- 
plitude, and the production of the pre-hadron occurs in 
the direct and conjugated amplitudes at different points. 
This also implies the appearance of an interference term 
between inside-outside production, which we consider in 
our calculation. 

The paper is organized as follows. In section II we in- 
troduce the kinematic variables and present the general 
structure of the amplitude of jet production. The ampli- 
tude is written in impact parameter presentation, which 
is especially convenient for the calculation of multiple in- 
teraction effects. 

In section III the cross section of inclusive hadron pro- 
duction is expressed via the light-cone Green functions 
describing the propagation of parton ensembles through 
a nuclear medium. The composition of these ensembles 
correlates with the coordinates of gluon-to-gg decays in 
the two amplitudes, direct and complex conjugated. Dif- 
ferent colorless parton ensembles propagating through 
the medium interact with different multi-parton cross 
sections. These cross sections are derived in Appendix 
B and expressed in terms of phenomenologically known 
cross sections of qq dipoles interacting with a proton. Us- 
ing known multi-parton cross sections one can solve the 
light-cone Schrodinger equations for the Green function 
describing propagation of partonic ensembles. These so- 
lutions are found in Appendix C. 

In section IV the cross section is presented as a sum of 
three terms, Eq. (32), the first two terms corresponding 
to the amplitudes for gluon decay both inside or out- 
side the nucleus, and the third term corresponding to 
inside-outside interference. These three terms are fur- 
ther evaluated employing different models for the pion 
wave function, which differ in the assumptions about the 



longitudinal and transverse momentum distributions of 
valence quarks in the pion. The unrealistic assumption 
of the Berger model that the quarks have no Fermi mo- 
tion results in a complete absence of nuclear effects, while 
more realistic models lead to a considerable nuclear sup- 
pression. 

Section V presents the main results and observations 
of this paper, as well as an outlook to future develop- 
ments. 



II. THE PROCESS -y'p nX IN IMPACT 
PARAMETER REPRESENTATION 



A. The amplitude 

In what follows we consider the reaction j*p ttX 
as an example of a hard reaction, neglecting the higher- 
twist terms calculated in [8]. Therefore, the space time 
development of this process and all the results are valid 
for any hard process producing a quark jet. 

The hard reaction 7*^ — > wX can be considered as a 
three-step process, as is illustrated in Fig. 1. In the first 
stage the incoming virtual photon knocks out a quark 
and a gluon of a bound nucleon, 7* + p ^ Qi + 5, which 
carry practically all the energy of the photon. We assume 
that Bjorken x is sufficiently large to neglect the contri- 
bution of the sea. It also allows to treat the gluon radia- 
tion process as incoherent. Strictly speaking one should 
integrate over the longitudinal coordinate of gluon radi- 
ation from —00 [12], however only a part of this path, 
Az ~ 1/xmjv, contributes coherently. We neglect Az as- 
suming that X is large. Thus, in the light-cone approach 
one can consider gluon radiation as instantaneous from 
the point of hard interaction, although the quark-gluon 
pair loose coherence at longer distances. One can come to 
the same conclusion analyzing the Feynman graphs. The 
corresponding space-time structure of DIS is studied in 
Appendix A. 

The second stage is the dissociation of the radiated 
gluon into a quark pair, 5 ^ 92 + 93, and the propagation 
of the colorless dipole 9291 through the nucleus. 

The third stage is the projection of the colorless dipole 
into a pion, 92 + 91 — > tt. The gluon and quarks prop- 
agating through the nucleus are assumed to experience 
only soft final state interactions (soft gluonic exchanges) 
with other nucleons-spectators, which usually cause at- 
tenuation. 

The amplitude of this process can be represented as. 
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Here $,r is the wave function of the qq Fock component 
of the produced pion; 



P2 = (1 - a)v-K - 



(3) 



p,r is the pion transverse momentum; k = (1 — a)pi — ap2 
is the relative transverse momentum of the q\ and qi in 
the pion; 



a = 



A = 



2 2 
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2E2 
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and £^1, E2, E3 and Eg are the energies of the three 
quarks and ghion. 

The longitudinal coordinates Zi are defined as follows. 
zi is the coordinate of the collision between the virtual 
photon and the nuclcon; Z2 is the longitudinal coordinate 
of the point of dissociation g 93^2; ^(j'i,'>^i^) is the 
amplitude of the process j*N qigXi, with the original 
impact parameters of the produced quark (f^f^ ) and gluon 
(rf ), at the point with coordinate Zi. 

The propagation functions (Green's functions) 
G{zf,Zin;ff,fin;E;{R}) in Eq. (2) describing the 
propagation of the fast quarks qi, q2, gs and the gluon 
in the medium, will be derived in Sect. III. Besides the 
initial {zimrin) and final {zf.ff) positions, they also 
depend on the coordinates {R} of the spectator nucleons 
with whom they interact via soft gluonic exchanges. 

The cross section of pion production off a nucleus is 
given by the amplitude squared and averaged over the 
positions of all nucleons in the nucleus, 



da 



d^p^d^p^dzh 



A I d%dzipA{b,z^) {\M\^)^^^. (5) 



Here Zh is the fraction of the photon light-cone momen- 
tum carried by the pion. We singled out the integration 
over the coordinates of the " active" nucleon participating 
in the hard collision with the virtual photon. 



effects become independent of Q^. Indeed, according to 
the uncertainty principle the smaller is the quark-gluon 
separation r, the faster they are expanding with trans- 
verse momentum kr ^ l/r, 



dr 
'dt 



1 

7e' 



(6) 



Here f = f^^ — f^; E = Ei + Eg is the total energy of 
the jet. 

If the initial size is small, then after a while its small- 
ness will be forgotten. Indeed, the solution of Eq. (6) 
reads, 



2/ X 2t 1 
At sufficiently long time intervals. 



AmNXBj 



(7) 



(8) 



where XBj = Q"^ /2miqE is the Bjorken scaling variable. 
At large the second term in (7) can be neglected, and 
therefore transverse size of the quark-gluon pair does not 
depend on any more. If this time interval is signifi- 
cantly shorter than the mean free path of partons in the 
medium, no dependence of nuclear effects should be 
expected. In cold nuclear matter the typical mean free 
path is several Fermi, so for t ^ 1 fm we expect a very 
weak dependence when xsj ^ 0.05. This condition 
is well satisfied in the region of Xsj ^0.1 dominated 
by valence quarks, in which we are focused. This effect 
probably explains the very weak dependence on Q"^ of 
nuclear ratios observed in the HERMES experiment [13]. 

In what follows we assume that xsj is sufficiently large 
to neglect the second term in (7), which is equivalent to 
the approximation in (2), 



r(ff,f«i) «r(ff) 5{rl 



r9i 



)• 



(9) 



III. GREEN FUNCTION FORMALISM FOR 
PROPAGATION OF PARTONS IN A MEDIUM 



B. <5'^-dependence 



At large photon virtuality the relative quark-gluon 



separation is small, 



l/Q, and the nuclear 



Although the hard reaction (DIS in Fig. 1) occurs on 
different nucleons incoherently, the multiple final state in- 
teractions of the produced partons proceed further coher- 
ently. Indeed, the mean transverse momentum squared 
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gained by a quark propagating through a heavy nucleus, 
as measured in the Dreh-Yan reaction at 800 GeV, is very 
smah of the order of Ap^ - 0.1 - 0.2 GeV^ [14], and is 
even several time smaller at lower energies [15, 16]. So 
a quark of energy Eg ~ 10 GeV interacts with coherence 
length Ic — 2Eq/Ap^ ^ 40 fm, which is quite long com- 
pared to the nuclear size. 

The cross section for the reaction 7*p — > nX corre- 
sponds to the product of the direct and conjugated am- 
plitudes, presented graphically in Fig. 2. The two am- 




OO 



1 



Ggiz2, zi;firf; E^; {R})G;{z,, zi; s|sf ; i^^; {R}). In 
the case of in-medium interacting partons, the Feynman 
path-integral representation is more appropriate. It al- 
lows to introduce differential equations for products of 
the Green functions and their conjugated functions, av- 
eraged over nucleon coordinates in the nucleus (see in 
more detail in the Appendix of [11]). Such a combina- 
tion, describing the propagation of the gluon and quark 
qi between points zi and Z2, has the form. 



Wi = {^Gg{z2,z,;rir,';Eg-{R}) 
X Gg, (Z2, ^i; ^ ; E,;{R}) G*g{z2,z,; s^sf- Eg- {R}) 



X G*g^{z2,zi-JTs1':EiAR}) 



(11) 



{R} 



which is a part of the final equation (18) (see below). It 
satisfies the equations. 
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2En 



■ PA {b, Z2) Si {f^^ , s| , ' , r| ^ 



dWi 
dzi 



Arf 

^e7 



■ pA{b, zi) Y.i{f{ 



Arf 

2e: 



Asf 

~2e: 



Asf 

2ir 



(12) 



(13) 



FIG. 2: Graphical representation for the direct and conju- 
gated amplitudes of the process 7*p irX. The DIS hard 
process occurs incoherently on the same nucleon with lon- 
gitudinal coordinate 21. The radiated gluon decays to g2?3 
coherently in the two amplitudes at the points Z2 and 23 re- 
spectively. The colorless pre-hadrons (dipoles) (/2?i and ^291 
created at Z2 and 2:3 respectively, are projected to the pion 
wave function in each of the two amplitudes. 

plitudes correspond to different impact parameters, (r^) 
and (si), of the participating partons, and different lon- 
gitudinal coordinates, Z2 and 2:3, of g ^ qq decay. 

One can see that specific partonic ensembles propagate 
through different intervals of the longitudinal coordinate. 
It is {2gqiqi} in the interval zi — Z2, {.^"Zs 92 9191} in the 
interval Z2 — 23, {91 91 92 92 93 93} in the interval Z3 — Z4. 

In the case of free propagation in vacuum the Green 
functions introduced in (2) have a simple form. 



G{z2,zi;f2,'ri;E) 



iE 



X exp 



27r(z2 - zi) 
iE{r2 - n f 



e(z2 - zi) 



2(2:2 - zi) 



(10) 



In equation (5) the amplitude (2) squared contains bi- 
linear combinations of Green functions with the same 
initial and final longitudinal coordinates, but differ- 
ent impact parameters. For instance the product 



with initial conditions. 



Wi 



Z2<Zi 

Wi\ = (5(f| 







ff ) 5(; 



?9i 



ff) 



(14) 



Here 6 is the impact parameter of the virtual pho- 
ton, and the nuclear density is normalized to one, 
/ cPbdz pA{b, z) = 1. El is the total cross section of a 
4-parton colorless system gqi gqi interacting with a nu- 
cleon target. It is important to notice that gqi and gqi 
are in color triplet and anti-triplet states respectively, 
while gg and qiqi are color singlets. 

The 4-body cross section Si can be represented as a 
linear superposition of elementary dipole cross sections 
of interaction of a colorless qq dipole with a nucleon, for 
which there exists a well developed phenomenology. A 
derivation presented in Appendix B results in the ex- 
pression. 



Ei(r'?SsS,s"?Sr^') = - 



+ agg{r'^^-s''^) + lagg{r^-s^) 



(15) 
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The solution of Eq. (12) with a rcahstic dipole cross 
section, and which is valid for large qq separations, can 
usually be obtained only muiierically (see [17]). However, 
in the reaction under consideration, with a highly virtual 
photon, the typical separations are small and one may 
rely on the approximation [18], 

(Tgg{n-r2) = Cin-f2f, (16) 

where the factor C is known from phenomenology. It 
depends on energy and should be properly chosen de- 
pending on the energy of the photon. With this ap- 
proximation and a constant nuclear density, pA{b,z) = 

I 



X 
X 

X 

X 
X 

X 



Here r and sj* are the transverse coordinates of the 
quark qi {i — 1, 2, 3) at the point with longitudinal coor- 
dinate Zj {j = 1, 2, 3, 4) for the direct and conjugated am- 
plitudes respectively (same for the gluons); a and a are 
the fractional momenta of the quark qi within the pion in 
the direct and conjugated amplitudes respectively. Cor- 
respondingly, the transverse momenta of the quarks qi 
and 92 are, 

Pi = apT^ + k; p[ = api^ + k'; 

P2 = (1 - a)p„ -K\ = (1 - a)pT, - k'] 

Ps = Ps; (19) 

The energies of the participating quarks and gluon in the 
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po 'd(R\ — 6^ — z^), the equations Eq. (12) can be solved 
analytically (see Appendix BB.l). 

To progress further, for propagation of the quark qi we 
employ the relation, 

G{z3,zi;f3,fi;E-{R}) = j d^r2G{z3,Z2;r3,f2;E-{R}) 
X G{z2,zi;f2,n;E;{R}), (17) 

which assumes that zi < Z2 < Z3. Then the cross section 
Eq. (5) can be represented as. 
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two amplitudes read. 
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= OiZhE; 




E2 


= (1 - a)zhE; 




Eg 


= {l-azh)E; 




El 


= azhE; 




E2 


= {l-a)zhE; 




Ej 


= {l-aZh)E; 






= E3 = {1- Zh)E. 


(20) 


The functions W2 


and W3 in (18) describing 


; the propa- 



gation of partonic ensembles through the intervals Z2 — Z3 
and Z3 — Z4 respectively, are defined similar to Eq. (11), 



dcr 



A 



X Z4^00 



hm 2 Re db dzidz2dz3 pA{b, Zi)Q{z3 — Z2)Q{z2 — Zi) e 



i{A-\-io)z2—i{^—io)z3 



d'^pzd'^p^dzh 
j d^rfd^rf'dhldhf r{ff ,f^') 

I d^ld'rfd^sld^sf Wi{z2, zv,rlfr ; si, sT ; rf , ; sf, sT ; Eg, Eg, El, El) 

J d-'rfd^rfd'rf d^s^dhf ^^2(^3, ^2; r§\rt > ; ^i, ; r|, ; s|, ; Ei, E2, E3, Eg, Ei) 
j d\f d\f d^rf dhf dhf d^sf 

W3{zi, Z3; ff , r-f , ff ; sf , sf , ; , rf= , ; 4, , ■,Ei,E2,E3,Ei,E2, E3) 



J d^Kd'^K'dada^^{a, k) (a, k') exp[ip[ ■ si' + ip2 ■ s^^ + ip3 ■ sl^ — ip[ ■ f^' — ip2 - 



«2 



ipk-rr\ (18) 



W2 = (Gg, [z3, Z2; ft , rr ; E3; {R}] Gg, [z3, Z2; rT , ; E2; {R}] Gg, [z3, Z2; ft , rt ; Er, {R}] 

X Gl^[z3,Z2-,4\^'-,Ei-,{R}]Gl[z3,Z2-,si,4-,Eg-{R}]\ (21) 

/ {R} 
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X [zi, zs; st , sT ; ^3; {R}] [z4, Z3; sf , ; ^2; {^}] G; [^4, ^3; s*!^ , si' ; ^1; {^}] 



(22) 



They are the solutions of the foUowing equations, 



with the boundary conditions, 
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2E. 



Arf 
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The function S2 is the total cross section of interaction 

of the colorless parton ensemble qi,q2, 93 Qi , 9 with a mi- 
cleon, where the pairs qiq2 and qi, qi arc each in colorless 
states, the pair <72.<73 is a color octet, and the pair qi,g 
is an anti-triplet. Correspondingly, S3 is the total cross 
section for the ensemble , 52 , 93 , 9i , 92 j 93 , where each 
pair qi,q2, 91,92, and 91,91 is colorless, while the pairs 
92,93 and 9293 are color octets. These cross sections are 
derived in Appendix B and have the form. 



52 = a,-g(f«= - f«i) +a,-g(r"'i - s «0 - agg(f«^ - s^^ + ^[Tqqif"' - f') - ag-g(f«^ - f«^) - agg(f«^ - s"')] 

o 

+ ? [(^Qgis"' - s') - (Tggir'^' - s^) + agg{f^' - s«) + a,-g(f - s»)] ; (29) 

o 

53 - ^,-,(r«^ - r"?^) + a,-,(s-«i - s«^) - a,-,(f«i - s«^) - a,-,(s"?i - r"?^) + a,-,(f«i - + a,-q(r'?^ - s"?^) + 
+ <Jqg (r - s-^^ ) + I [agg (f «^ - ) + a,-, (s «^ - s «^ ) + a,-, (r «^ -««=') + a^, (s «^ - f ) - agg (f - s"^ ) 

o 

- aggif"^' - s«^) - ag-g(f «^ " f ) " ^^^(s'^ - s ) ] . (30) 



I 

The equations (12)-(13) and (23)-(26) have been solved in Appendix C, in the approximation of Eq. (16) and for 
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a constant nuclear density. 



This amplitude can be represented as, 



IV. THE THREE PARTS OF THE CROSS 
SECTION 

In order to discriminate between production of the pion 
within or outside the nucleus, we rely on the approxima- 
tion of constant nuclear density, pA{b, z) = po'd{L^ — z^), 
where L = \J R"^ — b"^ . Then we split the amplitude equa- 
tion (2) in two parts, M = Mi + M2, corresponding to 
pion production outside, or inside the nucleus, i.e. Z2 - 
integration within intervals L < Z2 < Z3 and Zi < Z2 < L 
respectively. Correspondingly, the cross section equa- 
tion (18) contains three terms. 



Cr = (71 + (72 + (73, 



(31) 



which are the amplitudes Mi, M2 squared, and their in- 
terference, respectively. These terms correspond to the 
following splitting of the integrations over Z2 and Z3 in 
(18), 



Z4 Z4 Z4 Z4 Li Li L Z4 

j dz2 J dz3 = j dz2 Jdz3+Jdz2 J dz3+ J dz2 J 



dz-t 



(32) 

In what follows we consider cross sections integrated 
over transverse momenta of the pion and recoil quark, 



dui f 2 
-1—= dp,, 
dzh J 



d'ps 



ddi 



d^Pv d^P2, dzh ' 



(33) 



where i = 1, 2, 3. 

Later, the results of numerical calculations will show 

that the interference term is negative, (73 < 0. This 
can be understood on a much simplified example of an 
"empty" nucleus, i.e. free propagation of particles. In 
this case the amplitude of the fragmentation process 
I1 ~^ (I1Q2Q3 is proportional to the value. 
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Zl 



-i(A — io)z 



X Gg{z,zi;pg)Gg,{z+,z;pi)Gq^{z+,z;p2), (34) 
where =Pi +P2; 
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where 



M = Um e*^ / dz e 



m\^p\ ml+pl {pi+p2f 



(36) 
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2E„ 



2Ei ^ 2E2 + 



^ ^2 



(pi + P 2) 
2E„ 



■Zl. 



io; (37) 
(38) 



Now we can split the amplitude into two terms, M = 

Min + Mout, corresponding to gluon decay inside (zi < 
z < z) and outside {z < z < z+) the nucleus respectively. 
Then from (36) we get, 



Min = (e-**'' - e-"^^^) ; (39) 
Mont = ^ {e-'^'+ - e-'^') . (40) 



At 2;+ ^ 00 the first term in A4oui. vanishes, e *''^+ 
0, because of the imaginary term in i9, Eq. (37). So we 
get, 



1 

^' 

2 r 

2Re {MinM:^t) = -A|i_cos[i?(z-zi)]} 

= -|7W^n|^ 



(41) 

}; (42) 



(43) 



Thus, we conclude that the interference term (43) is nega- 
tive and exactly cancels the inside production term. The 
cross section in this case is given solely by the outside 

production. 

Of course, these simple results are valid only for 
hadronization in vacuum ("empty" nucleus). Presence 
of a medium breaks down these simple relations and 
makes the calculation of different terms in the cross sec- 
tion Eq. (31), performed below, much more complicated. 
Nevertheless, the negative sign of the interference term 
will be preserved. 



A. Pion production outside the nucleus 

We start with the first term ui, which dominates at 

high energy E, and is the easiest one to calculate, since 
in this case the functions M^i,2 contain just products of 
Green functions for free propagation of quarks in vac- 
uum. So the integration over longitudinal coordinates in 
(18) and transverse momenta in (33) can be performed 
analytically. 
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L 

ApoE^Fi'{x,Q^) J dzij dVf d^gs f(ff ) f(sf ) j S4(frfSsl<fsfdada 



dzhd?h ' ' - - - - J, 

X {L, zr, f|, r-f , si sr ; , ff, sf, sf )^(rf , ; s|, sT ), 
with the new notation, 



X (5 ( af2 — — — p — aS2 
1 - azh 



OtZh 



(44) 



= ^(^2 - ^2) / d^pd\Ko{mgp) KoirUgT) frs + p) K (^2 + A 

J \ 1- azh J \ 1 - azh J 



(45) 



^91 



Ri = 
Si = 



E, + Ex 



(46) 



(i = 1, 2) are the intrinsic separations in the qi— g pairs, 
and the coordinates of their centers of gravity, respec- 
tively; Eg+Ei=Eg + Ei= E. 

We also introduce the following combinations. 



Ra = Ri — Si. 



(47) 



The Jacobian for transition to the new coordinates is one, 

so. 



dV? dV? d^r^' 



d^Rt d^R^ d^r2d^S2. 



(48) 



It turns out that the function Wi factorizes in the new 
coordinates (see Appendix C), 

W, {L, zr, f|, ; s| , sT ; rl, ; sf , ) 

2 

F (^L, zi;R2,r2, S2; -Rf, n, si^ . 



E 



X exp 



2tt{L - zi) 
' _iE 
L — zi 



(49) 



Taking also into account that in (44) ff = ^ and sf = 
(i.e. ri = si = 0), we arrive at the relation. 



d^rf dhl d^R+ d^R2 f (ff ) f(sf ) 

X Wx{L,z,;?i,r'l\n.n':r!,r!,nin)^{R2 - S 
= j d^R+ |f p F{L, zr, r2, S2), (50) 



where 



F{L,zi\f 2,82) 

F (l, zi;R2, f2, S2; R^,fi,si^ ^- = ^ - = (^1) 



r 1 = s 1 = 



Thus the cross section equation (44) gets the form, 

L 

dai 



dzhd?b 



ANzl{l - Zhfpo j dzi j d^r2d^S2dada 



X F{L,zi\f2,S2;a,a)^{r2,S2), 



(52) 



where 



^{r2,S2 



d^pd'^T Koirugp) Ko{mqT) 



X $^ a; f2 + 



X 5 ( ar2 



Zh 

1 - azh 
l-a ^ 



p] a; S2 



1 - Z/l 



1 - azh 



p - as2 + — T 

1 - azh 1 - azh 



and 



N = 



f(f) 



(53) 



(54) 



In order to simplify the calculations we assume a 
factorized form of the pion light-cone wave function, 
<I>7r(a,r) = Lp{a)(f){r), and a Gaussian dependence on 
quark separation. 



(f){r) oc exp 



(55) 



where ^ is related to the mean pion charge radius squared, 

Perturbative fragmentation of quarks to pions in e^e 
annihilation and DIS was calculated by Berger [7] in the 
limit of (1 — z;j) <C 1. The pion wave function was max- 
imally simplified assuming that 95(0:) = S{a — 1/2) and 
fixing at zero the relative qq momentum. This simplifies 
the calculations considerably, since the function F can 
be obtained analytically, 



F{L, zi;f2,S2; a, a)a=a 
2 

exp 



2tt{L-Zi) 



I e 



2{L-zi) 



(ri - si) 



i po C{zh, a){f2 - S2f{L - zi) 



(56) 
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where e = azh{l — azh)E; 

C{zh,a) = c(l + a^zl 



4 



(57) 



Notice that in this case the expressions for W2 and W3, 
equations (21) and (22), also are much simphfied. 

We can perform the integration over the transverse co- 
ordinates and momenta using the integral representation 
for the modified Bessel functions, 



We use here the following notation, 

Di = {l — aft'w + m^{l — azh)'^u; 

D2 = (1 — Q:)^t W + TO^(1 — Q;2/i)^m; 

D3 = COS [uji{L - zi)] COS [lo2{L - zi)] ] (62) 



2 — eeu]iui2 cot [wi(L — zi)] cot [uj2{L — Zi)] ; (63) 



i 00 

Ko{mqp) Ko{mqT) = J 



dw 



X exp 



2w 



-1 
^2 



1 + V 1 — V 



w 



(58) 



Then in the case of equal sharing of longitudinal momen- 
tum by the pion quarks we arrive at a simple result. 



d(Ji 
dzhd'^b 



Li 00 

AN zl{l- Zhfpo j dzi J dw 



2(1 - ^1^2) 

- SWij'^ COt[L0i{L - Zi)]} , 



{£UJ20^ cot [uii{L - Zi)] 



(64) 



where (3 = a — va; 7 = a — jia; e = Eazhil — cczh); 
£ = Eazh{l - azh); toi = \/-iXi; 102 = V-tM'-, 



Ai,2 — C po 



2e£ 



sa + sb— \/ [ea + £6)^ — ^eec?- , , 
= £v = ~~2c ' ^ ^ 



dv 



rrig + aD{w,v, zi) ' 



(59) 



where 



D{w,v,Zi) = {l — v^)w 



+ {1 — x) m \ ua 




+ 2iv{l-Zh) 



mi{L-zi) 



(60) 



Notice that although the function D{w,v,zi) is com- 
plex, the expression (59) is real. 

In the limit of ^ ^ in Eq. (55) the cross section 
Eq. (59) does not depend any more on the interaction 
with the medium, which is characterized by the con- 
stant C. Thus in the Berger model for fragmentation 
{a = 1/2; (k^) = ^ = 0), the interaction of the quark 
and gluon with the medium does not affect the value of 
the cross section (Ti,Eq. (59), and only modifies the trans- 
verse momentum distribution, which is an effect beyond 
the scope of this study. 

In another hmiting case $77(0;, k) = ^{a) 5{k) the cross 
section gets the form, 



da I 
dzhd'^b 



= AN zf^{l — Zhf Po j dzi J dw J da J da 



-LOO 

(1 - azh){l - azh)'fi{a)ip{a) 



Di D2 D3 



(61) 



a = (1 - azhf + - azh; 

9 

b = (l-azh) +-azh; 

9 

c = {I- azh){l- azh) + -{azh + azh). (66) 

Nuclear effects for this part of the cross section, cti, are 
shown in Fig. 3 in the form of ratio, 



Ri{zh) = 



dai/dzh 



dai{C = 0)/dzh' 



(67) 



where both the numerator and denominator are the cross 
sections on the nucleus integrated over impact parameter, 
however in the denominator we eliminate the influence of 
the medium fixing the imaginary part of the light-cone 
potential C = 0, so the quark and gluon propagate like in 
vacuum. The nuclear cross section, here in the numerator 
and in what follows, is calculated with C = 3. This value 
agrees with extrapolation of the saturated cross section 
[12] down to medium high energies, as well as agrees with 
data on nuclear broadening of transverse momentum [19] . 

We performed calculations for four cases: 

I. ^jrict, k) fx S{a — 1/2) 6(k) (Berger approximation); 

II. $7r(Q;, k) = (p{a) (5(k), where ipia) oc a{l — a); 

III. $^(a, k) = 6{a — 1/2) (/)(k), where </)(k) oc 
exp(-^KV2); 

IV. ^Tria, n) = ip{a) 4>{k), where <f{a) oc a(l— a) [20, 21], 
(/)(k) cx exp(— ^k^/2). 

Comparing curves I, II with III, IV in Fig. 3 one 
can conclude that the transverse motion of quarks (^ = 
(k^) ^ 0) significantly affects the nucleus-to-proton ra- 
tio. At the same time, the a-distribution, i.e. longitudi- 
nal motion of quarks in the pion, has almost no influence 
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N 




^ 0.9 



FIG. 3: Comparison of ratios Ri{zh), Eq. (67), for pre-hadron 
production outside the nucleus, for different models for the 
pion light-cone wave function. The four variants I-IV (see 
text) differ by absence or presence of longitudinal and trans- 
verse motion of the valence quarks in the pion. The calcula- 
tions are done for lead at B = 10 GeV. Since the recoil quark 
should be ultra-rolativistic in order to rely on the Green func- 
tion method, we restricted the range of Zh < 0.9. 



on the nuclear effects. Indeed, the curves III and IV 
arc nearly very close to each other. In what follows wc 
assume that a2,3 are also insensitive to the form of the 
a-distribution, so we will continue our calculations in the 
approximation III, ^-jria, k) = 6{a — 1/2) exp(— ^k;^/2). 



Here 

if2 



- /9o[(ci - C2)(-22 - + (C3 - C'2){zs - Z\f 

+ C3{L- Zif]+2ie\{2x + \){Zi- Z2) 

- 2iSw{ [{x + Xf{z3 - Zif + X^{Z2 - Zif] 

X cot [w (2:3 - Z2)\ 

2{X + \)\{Z3 - Zi){Z2 - Zi) 

sin[w(z3 - Z2)] 
+ m-xf{L-Zr)\ 



and 



(69) 



X 

X = 

P = 
Q 

I = 

C2 = 

C3 



, xE 

9. 

P 

4-17X + 22a;2 
4(l-a;)2 ' 

4 - 8.T - 5a;2 



4(1 -x) 
^, X 
.1+4+"- 



'■) C; 



4 p 



4(1 - xfC. 



(70) 



One can see that even in the limit ^ ^ the cross section 
Eq. (68) is still sensitive to the constant C due to the 
presence of sin[u;(z3 — Z2)] in (69). 
The result for the ratio 



R2{Zh) 



d(j2/dzh 



da2{C = Q)/dzh 



(71) 



B. Pion production inside the nucleus 



calculated for lead at 10 GeV is depicted in Fig. 4. As 
one could expect, this contribution is more suppressed, 
since the gluon decays inside the nucleus producing the 
colorless pre-hadron which propagates and attenuates in 
the medium. 



The second term in the cross section Eq. (72), which 
corresponds to gluon decay inside the nucleus, has the 
form, 



da2 



dzhd'^b 



1^ 1^ 1^ 
irANpoKe J dzi j dz2 j dz^ 



-L Zl Z2 



—i£(jexp —^^{z3 — Z2) 



E^(l + ^ H2) sm[u{z3 - Z2)] ' 



(68) 



C. Interference of amplitudes JVfi and M2 

Eventually, the third term in Eq. (31), corresponding 

to interference of the two amplitudes, after integration 
over transverse variables gets the form, 

L L 00 

-L Zl 



(72) 



—i£uj exp 


-^(^-^2) 








sin[u;(L — Z2)] 
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FIG. 4: Lead-to-proton target ratios for the cross sections 
corresponding to pre-hadron production outside {R\) or inside 
{Jh) the nucleus in both ampUtudes. The ratio for inside- 
outside interference is shown by the curve indicated by Rz- 

The ratio of the full cross sections ai + a2 + 0-3, Eq. (76), is 
depicted by the solid curve indicated by R. Notice that the 
interference term gz is negative. 



Here 



Hi 



m; 



+ 4^(l-x)2(L-zi)2 + 2i5A'(L-Z2) 
- 2i£:w| [(L - zif + (z2 - zxf\ cot[a;(L - Zi)\ 
(L - zi)(z2 - zi) ' 



sin[w(L — z-i) 
,2 



2ie{\-Zh){L-zx)\ (73) 



m„ 



w 



■ 2i£L0 cot[w(L — Z2)]; 



(74) 



B3 = !!^(l-:,.)(i-^l) + 2^fV-^ 
ws \ 1 — a; 

— 2i£ui\^L — z\) coi[ui{L — Z2)] 



-22 - Zi 



sin[w(L - 22)] J ■ ^^^^ 

The values of R3{zh), defined similarly to Eqs. (67), 
(71), are depicted in Fig. 4. The solid curve in this figure 
presents the final results for the ratio of all terms in (72), 



R{zh) 





1 da-2 


, das 


\ dzh 


dzh 


dzh 


dai _|_ 


da-2 1 


da:i \ 


dzh 


dzh 


dzh J ( 



(76) 



Fig. 3 docs not contain information about relative con- 
tribution of different terms in (76) to the cross section. 
To show that wc depicted the fractions di/iai +a2+ tra) 
in Fig. 5. All three terms in the numerator of (76) are 




FIG. 5: Fractions ai/{a\ + 02 + crs) calculated for lead at 
E = W GeV. 

of the same order, but the last one presenting interfer- 
ence, is negative. The latter was expected according to 
the calculation performed in Sect. IV for hadronization 
in vacuum. We see that at 10 GcV the fractions of the 
cross section corresponding to production inside and out- 
side the nucleus are about equal, while the former is more 
suppressed according to Fig. 4. This is, however, a clas- 
sical interpretation, the inside-outside interference term 
0-3 does not allow to classify events this way. 

The nuclear effects represented by the ratio R{zh) de- 
pend on the photon energy, and the higher the energy is, 
the weaker is the nuclear suppression. This is the obvi- 
ous manifestation of color transparency [18]: the initially 
small quark-gluon separation (see Sect. TIB) is evolving 
slower at high energy due to Lorentz time dilation. The 
energy dependence is illustrated in Fig. 6 by some exam- 
ples. 



V. MAIN RESULTS AND OBSERVATIONS 

This paper presents the first attempt to describe 
hadronization of a parton propagating through a medium 
on a fully quantum-mechanical basis. For hadroniza- 
tion in vacuum we employ the Bcrgcr model [7] of per- 
turbative fragmentation, improved in [8]. This mecha- 
nism, imbedded in a nuclear environment, is illustrated 
in Fig. 1, and the associated space-time development in 
Fig. 4. We employed the path-integral formulation [22] 
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FIG. 6: Ratio R(zh) defined in (76) for lead at different pho- 
ton energies E = 5, 10, 15 GeV. 

of quantum mechanics, which describes propagation of 
partons and partonic ensembles in terms of the light- 
cone Green function formalism. This technique properly 
includes all phases and takes care of all coherence phe- 
nomena, including formation of hadronic wave functions 
and color transparency. 

The important observations of the paper can be sum- 
marized as follows. 

• Contrary to the usual expectation based on clas- 
sical intuition, even if the radiated gluon always 
decays outside the medium and the produced ^293 
pre-hadron has no final state interactions, there is 

a considerable nuclear suppression for pion produc- 
tion, as is demonstrated by the ratio Ri{zh) in 
Fig. 4. Notice that no energy loss effect or final 
state absorption arc involved in this result. The 
suppression is caused by multiple interactions of 
the partons in the medium affecting the overlap of 
the pre-hadron and pion wave functions, even if the 
pre-hadron is produced far away from the nucleus. 
The effect of nuclear suppression is subject to color 
transparency and is controlled by the size of the ef- 
fective dipoles. The latter is evolving starting from 
a very small separation ^ l/Q in the hard reaction 
initiating the jet. The magnitude and energy de- 
pendence of nuclear suppression is similar to what 
is known for clcctroproduction of p-mesons in the 
regime of short coherence length [23, 24]. In that 
case a dipole is also produced in a small-size con- 
figuration and then is evolving with a speed depen- 
dent on energy. 

• Although much more involved, the effect of nu- 



clear suppression of prc-hadrons produced outside 
the nucleus is in some respect analogous to gluon 
shadowing (there are also important differences). 
Indeed, in the nuclear rest frame gluon shadowing 
looks like suppression of gluon radiation by multi- 
ple interactions [11, 12] (Landau-Pomeranchuk ef- 
fect). In this case there are no colorless objects to 
be absorbed in the nucleus, yet the production rate 
of gluons is affected by the medium. In this case 
gluons radiated inside and outside the nucleus also 
interfere. 



• The novel feature related to the quantum- 
mechanical treatment of the problem, is the pro- 
duction of the pre-hadron both inside and out- 
side the nucleus. This is analogous to the Twin 
Slit Interference Experiment in quantum mechanics 
when a particle propagates simultaneously through 
both slits. Interference of the amplitudes with in- 
side/outside pre-hadron production has a consider- 
able effect on the nuclear absorption. This inter- 
ference term in the cross section is large and nega- 
tive, as is explained in Sect. IV on the example of 
hadronization in vacuum. It is not a surprise that 
the possibility of pre-hadron production inside the 
nucleus leads to more suppression due to attenua- 
tion of the colorless pre-hadron g2?i- 



• Suppression of hadrons should be much stronger 

in the case of a dense medium created in heavy 
ion collisions. This effect is completely missed in 
calculations based on the energy loss scenario [25]. 
In fact, it should account for a substantial part of 
high-pT hadron suppression observed in heavy ion 
collisions [26]. This may also explain why the ob- 
served suppression, when is related solely to energy 
loss, demands an unrealistically high density of glu- 
ons radiated in heavy ion collisions [27]. 

While the performed analysis highlights the novel fea- 
tures of in-medium fragmentation brought by a rigorous 
quantum-mechanical treatment of the process, it is still 
not sufficiently realistic to be compared with data. Frag- 
mentation was calculated in the Born approximation, and 
the main lacking element is vacuum energy loss due to 
gluon radiation caused by the initial hard interaction [28] . 
Such a modification is expected to shrink the distances 
Z2 — Zi and z^ — Zi and make them Q^-dependent. More- 
over, vacuum energy loss caused by gluon radiation leads 
to a distance for pre-hadron production which vanishes 
in the limit Zh ^ 1 [2, 4, 29]. Energy conservation also 
causes nuclear suppression toward the kinematical limit 
Zh = [30]. These corrections may only enhance the 
statements listed above. We plan to work on this prob- 
lem and publish elsewhere. 
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Appendix A. SPACE-TIME STRUCTURE OF 
THE DIS VERTEX 

The ampUtude of the process Iqo ^''Zi92<Z3 has the 
form, 



M = 



■(I) jW 

' 



(A.l) 



where j^^ and jj!^^ are leptonic and hadronic currents 
respectively. The latter can be presented as a sum of two 
terms, 



j{h) _ j{a) , j{b) 



(A.2) 



corresponding to graphs Fig. 7a and b, respectively. The 

y. / - y* 





a b 

FIG. 7; Two Feynman graphs for reaction 7*gi q'lq^qa- 
amplitudes have the form, 



+ P2 +P3h^Lu{po) 



« -^u{pihaG{pi 

X "(^3)70^(^2); 

X u{p3)'JaV(p2), 

where = (^2+^3)^- 

For a massless quark, mq = 0, 



+ 10 



(A.3) 
(A.4) 

(A.5) 



Correspondingly, 



(A.6) 



The propagator G{x2 — xi) describes propagation of a 
quark through the interval X2 — xi absorbing the virtual 



photon at one point, and radiating a virtual gluon at 
another point. 

One can obtain an amplitude of the reaction Iq^ — s- 
I'lrq^ from the amplitude Eq. (A.l) projecting the pro- 
duced g29i pair to the pion wave function. Generally, 
this is quite a nontrivial problem, which, however, much 
simplifies in some approximations. In the Berger model 
of a " frozen" pion [7] one neglects the intrinsic motion of 
the quarks in the pion, p2 — pi — \pTn which is certainly 
not a realistic approximation. In this case the projection 
is fulfilled using the relation [8], 

^u(^)(pi)t-(-^)(p2)sgn(A) = \l5{Pn + ra^). (A.7) 



Another approximation is to neglect the quark and 
pion masses, = m-r; — 0. In this case, 

^zl(^)(pi)f(-^nP2)sgn(A) = Va(l"a)75P.r, (A.8) 



where a is the fractional light-cone momentum of one of 
the valence quarks in the pion. Then the projections of 
the components of the hadronic current to the final Trqa 
state have the form. 



= / da -^^^fer^u(p3) 75^776(^3 +P7r) 7m" (Po); 



M2 



(A.9) 



J^!"^ ^ I da ^fri'^'^ u(p3)75G'(ap^ - Q)^ ^,p.^u{pQ) , 



(A.IO) 



where = [pa -I- (1 — Q;)p^]^. 

Since according to the Dirac equation, ?2(p3)p3 = 0, we 
can replace the product p-nG{pT^ +P3) in the component 
Eq. (A.9), by 

G(p^ +P3) = (P^ +P3)G(p^ +P3) = 1. (A.U) 
Then the effective propagator. 



G{X2 - xi) = 



1 



{2n) 



d^pGip) e*P(^=^-^i) =5(x2-2:i). 



(A.12) 

is not zero only when X2 — xi, i.e. the virtual gluon is 
radiated at the same point where the virtual photon is 
absorbed. 

For the second component of the hadronic current, 
Ji^\ Eq. (A.IO), it was demonstrated in [8] that 
its longitudinal-to-transverse contribution ratio is sup- 
pressed as aL/cTT ^ pjj_/(l ~ z)Q^- Therefore, we will 
keep only the transverse part of the current J^^ . 

Then, using the kinematic relation, p^^ — z^ps/il — 
Zh) + 0{pTr±) and applying the Dirac equation, u(j>3)p3 = 



14 



0, the transverse part of the hadronic current Eq. (A. 10) 

can be modified as, 



u{p3)'y5G{apT, - Q)^fj,p„u{po) = G u{p3)^5l±u{po) 

Q 



+ O ( 4f ) ■ (A-13) 



where 
G = 



Zh 



Q'^Zh{l - azh) - ap^j_ I - azh 



O 



r 

(A.14) 

Thus, neglecting corrections of the order of {p^^/Q"^ we 
arrive at the relation between the two components of the 
hadronic current, 



'•^ - I- azH ■ 



(A.15) 



Therefore, the effective quark propagator for the sec- 
ond component of the hadronic current, corresponding to 
Fig. 7, should have the same property Eq. (A. 12) shrink- 
ing the interval X2 — x\. 



Appendix B. MULTI-PARTON CROSS 
SECTIONS Si 

The effective cross sections , which are linear combi- 
nations of dipole qq cross sections, can be derived within 
the Born approximation. The derivation is quite lengthy 
and not easy. It is much easier to used a set of equa- 
tion which correspond to different limiting configurations 
within the multi-parton state. The way how it works is 
explained further in concrete examples. 



B.l. 4-body cross section Si(f s*^^, f s") 

To proceed, let us start with Ei , which is the total cross 
section of a 4-parton colorless system gqi gqi , in which 
the gqi and gqi pairs are in color triplet and anti-triplet 
states respectively, while the pairs gg and qiqi are color 
singlets. Apparently, the system is symmetric relative to 
interchanges r ^ ^ and f'^ ^ s'^, so we get, 

+ bi [a^,{r9 _ gi) + ^.^^g9 _ ^--9)] + ciaq-,(f « - s«) 

+ diaggifs - sO). (B.l) 

This expression, at = f, must turn into the known 
3-body {g, q, q) cross section. 



leading to the relations, ai = bi + di = ^ and hi + Ci = 
— i. Otherwise, fixing r'^ = s'^, we expect. 



Hi^aggir^ 



-a,-,(f«-5«). (B.3) 



This leads to new relations, ai -|- 6i = and di = | . 

Thus, we arrive at the coefficients in (B.l), oi = —hi = 
|, ci = 1 and di = |, which proves Eq. (15). 



B.2. Cross sections S2(r f"'^, f , s s ''^) and 
S3 (f « , r « , r «i , s"^ , s " , s «i ) 

Again, relying on the symmetry relative the replace- 
ment f^swe can write. 



S3 = 


as [o-gq(r^' 






s«^)] 


+ 


63K-,(r-'«^ 


-f"^) + 


(^ni^"' - 


s«^)] 


+ 




-r«^) + 




S«3)] 


+ 


d3 K-,(f«^ 


-s^') + 




^92)] 


+ 


e3K-,(f«^ 


-s'^') + 






+ 


/3[^5,(r«"^ 


-s«^) + 


^mis""' - 


r«3)] 


+ 






-s«^) 


+ 











(B.4) 

We can simplify this expression by considering known 



limiting combinations. For r 



"3 and s «2 



g<i3 = gg we have to get E3 Ei(f9,f«i,sS,s«i). This 
condition leads to the following relations, 03-1-63 = —d^ — 
63 = |, 53 = -1, 2/3 + /13 + fej^ |. ^ 

The next possibility is to fix s = s = and r = 
s ^1 . Then we should arrive at the 3-body case, 

S3 ^ C73(f «^ , , s^) = \ [agg{sa - f 9^ ) 



)-^<J,q{r'''-r'''). (B.5) 



This results in additional relations, fs + hz = fz + k^ = ^, 
C3 = — I and as + = 63 + C3 = 0. 

Eventually, after fixing the coordinates differently, 
gq2 = ^93 = gg and f«i = r^^, Eq. (B.4) simplifies to. 



S3^(73(s«^f«^s«). 



(B.6) 



Correspondingly, new relations emerge, 03 + 63 = + 
h = §7 63 -I- C3 = ^3 + 53 = /J3 + /s + 4 + 63 = 0, and 
e3 = -i. 

Solving these sets of relations we get the coefficients in 
(B.4), 



^3 = -ds = 53 = /I3 = ^3 = 1, 

63 = -63 = -C3 = /3 = 77' 



(B.7) 



(B.2) 



which lead to Eq. (30). 

Eventually, one can get Eq. (29) for 

S2 (f f«^f 91,8 9,591) by fixing s^^ = s« = 
in Eq. (30) for S3 (r 9^ , f 9^ , f 9i , s 9^ , s 92 jqiy 
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Appendix C. FUNCTIONS Wi 
C.l. Pull calculation 

As we already mentioned, equations (12)-(13) and 
(23)-(26) can be solved analytically, provided that the 
nuclear density is constant, pAib, z) = po, and the dipole 
cross section has the simple form = Cr^. In this 
case the equations, which are bilinear in the interaction 
potential, can be solved following Rcf . [22] . 

We demonstrate here the method for the example of 
Eq. (12), which can be represented as, 



I— — =HiWi. 

OZo 



(C.l) 



The effective Hamiltonian Hi can be written as a sum of 
the effective kinetic and potential energies, 

Fi=Ti+yi. (C.2) 

T hen Wi can be presented as, 

Wi{z2,zi)=Ni{z2,zi)e'^\ (C.3) 

where 



Z2 

Si = j dzLi{z). 



(C.4) 



(C.5) 



The effective Lagrangian has the form, 
Li=Ti{z)-Vi{z), 
where the potential term reads 

Vi{z) = -ipoSi(fs(^),f«H^),s9(^),s«H^)). (C.6) 
The kinetic term has the form, 



Ti{z) = ^E,{v^{z)f + ^Ei{v''^{z)f 

- lE,{u^{z)y-lEi{u'^^{z))\ (C.7) 



where 



iJ9{z) = -^f\zy, 

dz 

iT9i(z) = -^r'^'iz); 
dz 

qi^{z) = -^fi^{z). 
dz 



(C.8) 



lutions of the Euler-Lagrange differential equations, 

d f dLi\ dLi 

= 0; 



dz \dv3 J df9 
d ( dLi \ dLi 



dz J dfi^ 

d f dLi \ dLi ^ 



0; 



dz \du3 J ds9 
d ( dLi \ dLi 



dz \ dull J as 91 



(C.9) 



It is convenient to use the coordinates of the center of 
mass, 



R 



E„r9 + Eifi^ 



Eg+El 

where x = azh, and 



= (1 - a;)r » + a;f«S (C.IO) 



S = = (1 - x)s^ + 



E„ + Ei 



(C.ll) 



where x = azh- 

The relative separations arc given by. 



The corresponding velocities read. 

Viz) = iRiz); 

IJiz) = ±Siz); 

v{z) = ^r(.); 

u{z) = ^^s{z). 



(C.12) 



(CIS) 



In the new variables the Lagrangian Eq. (C.5) gets the 
form, 

L,{z)^\E{y'-U'^)+\ev^ 

- ^e-u2 + lpoEi(f,s,i?-5), (C.14) 

where E = Eg + Ei = Eg + Ei = E. 

Then, we make the following combinations of the cen- 
ters of gravity coordinates, 

1 



and velocities. 



The transverse separations as functions of z are the so- 



R+ = ^-{R + S); 
R = R — S', 



V+ = ^(V + U); 
V- = V-U: 



(C.15) 



(C.16) 
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Notice that the cross section Ei Eq. (15), which en- 
ters the potential term of the Lagrangian, Eq. (C.6), is 
independent of Therefore, the Euler-Lagrange equa- 
tions (C.9) written via new variables R~^,R~, f, s have 
a simple solution, j£V~ =0, i.e. 



V-(z) = Const 

R-{z) = R~ + {z- zi)V-. 



(C.17) 



Then, for the first term in the Lagrangian Eq. (C.14) 
the integral Eq. (C.4) can be calculated as. 



E 



E 



Z2 - Zl 



j dzV + 

Zl 

{R+ - R+){R2- - Rn- 



V- 



(C.18) 



Thus, for this part of the integral Eq. (C.4) we did not 
need to know the explicit form of R~^{z), which is rather 
complicated. 

In order to calculate the rest of the integral (C.4), we 
need to know f{z) and s{z). The potential Eq. (C.6) can 
be represented as. 



Vi = 



af'^ - 2b f- s + cs'^ + 2df- R~ 

^ \ 2" 



- 2es- R~ + 



f{R-y 



(C.19) 



where 



= —iCpo [l + a^zf^ + azf^A'j 



1 + aazl^ + {a + a) 



-iCpo 



a 
b 
c 

d = —iCpo 



Zh 



iCpo [l + a^zl + azh^ 
1" 



/ = 



-iCpo 
-iCpo 



azh 



1 

azh + - 



(C.20) 



Then the Euler-Lagrange equations lead to the follow- 
ing linear equations for r[z) and s{z), 



f = —af+bs — dR~ 



s = —bf+cs — eR~ 



(C.21) 



where e = x{l — x)E, e = x{l — x)E. To make these 
equations homogeneous we switch to new variables, 

r' = r + vR~; 

s' = s + CR~, (C.22) 
where v and C are solutions of the algebraic equations, 

va — (b = d 

vb-Cc = e. (C.23) 
Then, f'{z) and s'{z) satisfy the homogeneous equations. 



e , „ = —ar 



dz^ 
d^s' 
dz'^ 



bs'; 

= -bf' + cs'. 



(C.24) 



The solution of these equation is, 

f'{z) = A sin(ti>iz) + B cos(ti>iz) 

-I- ^[C sm{LL!2z) + D cofi{uj2z)\] 

s'{z) = p[A sin(a;iz) + B cos(ci;iz)] 

-I- C sin(w2-2) + D cos{u!2z). 



(C.25) 
(C.26) 



Here 0^1,2 = \/M,2, where A1.2 are the solutions of the 
quadratic equation (a — eA)(c -|- eA) — 6^ = 0, and p = 
(o- eAi)/6, ly = {c + eX2)/b. 

The vectors A, B, C, D arc fixed by the boundary 
conditions r'(-^i,2) = ^1' 2 ^^'^ ^'(-^1,2) = ^1 21 



A = 

B = 

C = 

i) = 





1 pici2 


- P2C11 


1 


- pV S11C12 


- S12C11 




1 P\S\2 


- P2S11 


1 


- pv C\\S\2 


- C12S11 




1 nC22 


- T2C21 


1 


- pv S21C22 


- S22C21 




1 T1S22 


- T2S21 


1 


- pv C21S22 


- C22S21 



(C.27) 



where 



Tl,2{z) = Sl,2(^) - i^?^l',2(^); 



(C.28) 



and Sij = sin(ajiZj), Cj.j = cos(a;jZj), i, j = 1,2. 

Now we are in a position to perform the rest of inte- 
gration in Eq. (C.4), and we arrive at the final expression 
for the action. 
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Si — — — (^2''' ~ -^l'')(-^2 ~ -^1 ) + ^ ^ 



AZ12 

a;2-E'x(l — X 



2(1 - iiiy) 



(Pi +P2)cot(a;iA2;i2) 



2pi ■ P2 
sin(t<;i Az) 



+ 



2(1 - nu) 
Ex{l-x)C 

^Z12 



Ti + T2 ) cot(a;2Azi2) 
(s2 - si) {^R^ - ^1 



sin(a;2A2;i2) 

(-^2~ ~ -^r)^ 

2A2;i2 



-Ex(l — x)v 



(f2 - fl) - 



-(/ + ^;2a + c2c-2t;C&) 



(4-) + (4-) 



Ai;i2) 

£;3;(l-x)?;2+£;x(l-x)C^ 
^ \ 2 



A212. 



(C.29) 



where Azi2 = Z2 — zi] x = azh] x = oiZh- 
Now we can solve equation (C.l), 

i-^-Hi^ e'^' = i^Azu) e'^' , (C.30) 
where 



Besides, for the parameters defined in (C.21) e = e, and 
the parameters ^ 1. It turns out that it is more 
comphcatcd to perform a transition in the found solu- 
tion for the action Si{a — a 0), than to repeat the 
derivation specifically in this limit. 

In this case the Lagrangian (C.5) gets a simple form, 



^{AZ12) 



+ ui cot(a;iA2;i2) + u)2 cot(a;2Azi2). 



(C.31) 

Then, from (12) and the boundary condition (14) we find 
the factor in (C.3), 



Ni{Azi2) 



J. 



E 



E'^x{1 — x)x{l — a;)a;ia;2 



\ 27rA2;i2 / (27r)2 sin(a;i A2;i2) sin(a;2 A2:i2) 

(C.32) 

The derivation of the functions 1^2,3(22, -Zi) is analo- 
gous, but rather cumbersome, so we skip it here. 



C.2. Approximations 

1. Function W\ 

The expressions for Wi{z2,Z\) {i = 1, 2, 3) simplify, 
if X = a; (i.e. a = a), and the parameters a,... / in 
Eq. (C.19), which are functions of x and x, are related if 

a(x = x) = h{x = x) = c{x = x) 

d{x = x) = e{x = x) (C.33) 



Li{a = a) =EV+ -V- +Ex{l-x)v+ -v- 
\[a{r-f + f{R-f + 2ef-R-] , (C.34) 



where — ■§^r^, = ir + s)/2, f = r — s. 

From the Euler-Lagrange equations of motion it follows 
that -^V- = -^v - =0, so 



V- = 



V = 



R2 ~ ^1 

AZ12 ' 
AZ12 ' 



(C.35) 



and 



R = Ri +v {z — zi); 
f~ = + v~{z — Zi). 



(C.36) 



This is sufficient for calculating the action Eq. (C.4), 
and we arrive at, 



Si = 



E 

AZ12 
fAZi2 



[r+ - R+) (^2- - Rr) + x{l - x) {f+ - r+) (rT - rV) 

(i?2-)' + i?2"i?r + (i?r)^ 



aA2;i2 



(f2 )Vf2 fl + (fl )^ 



eAzi2 



2i?2 ^2 + 2-Rrfr + R2'r'{' + Rir^ 
I 



(C.37) 



Then the coefficient A''i(A2i2) in Eq. (C.3) gets the very simple form, 

Ari(A2i2)=a;(l 



x)( 



E 



V 27rA2;i2 



(C.38) 
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2. Function W3 

The next case is W3, which is simple due to the symme- 
try relative to interchange r s ^* . First we introduce 
the Jacoby coordinates, 



R = 




f = 






P = 




S = 




s = 


X2S^^ + 2:35'^ 


X2 + X3 


T = 


gQ3 _ gq2 



(C.39) 



Here xi = az,,, X2 ^ {1 - a)zh, X3 = 1 - Zh- 

We also introduce combinations of the Jacoby coordi- 
nates, R+ = \{R + S), R- = R - S, r+ = ^{f+s), 



The Lagrangian L3 can be represented as, 



r = r - s, = ^{p + f), p = p- f. 



L3 = EV+V- +x{l-x)E{)+v- 



where y± = ^J?±, = ^f±, a;± = 



a = — 

b -- 

c = 

d -- 



iCpo 

iCpo 

iCpo; 

iCpo 



1 , 22, azh- 



1 



9 (1 - a)zh{l - Zh) 
4 l-Uh 



, 1 9 {1 - Zh){l + azh) 



1 - azh 



= -iCpi 



azh + ■ 



(C.41) 



Again, the equations of motion lead to the relations, 

V- = —v- = —v-=0; 

oz oz 

Z - Z3 



dz 



1 



a r 



R- = 



P = 



i?j + {Ri — R2 
r,- + {R^-f2~) 
Pi + {Pi - P2) 



AZ34 
Z- Z3 _ 

2-2:3 



A2; 



34 



2^^v -y+b{p-f+c(R~ 

+ 2df~p~ + 2ef~R~ + 2fp~R~^, (C.40) These relations lead to the following action, 
I 



(C.42) 



24 

S3 = J dzLs 

Z3 



1 



AZ34 

AZ34^ 



)^ + +f2 r*! + (fi )^ + & (P2 y + +P2 Pi + [pl ) 



R. 



y ++R2-Rr +{Riy 



+ d [2f2 p2 + 2f-^ P2 +f2 Pl + ^2 P2 ] + e 2f2 R2 + 2ri R2 +f2 R^ + f2 ^2 



I 



2p2 i?2 + 2/?! i?2 + P2 ^1 + P2 ^2 



(C.43) 



Eventually, we arrive at, 

W3{z)=N3{Az34)e'^', 

where 

N3{Az34) = a{l -a)zl{l-Zh) 



(C.44) 



E 



27rA234 



(C.45) 



5. Function W2 

The calculation of is more involved. First we switch 
to new variables. 



R = 



r = 



P 
S 
s 



Xif^^ + X2r'^^ + X3f'^^; 
X2r"''- + -r-.ir'''-' 
X2 + X3 

rl3 _ T^'Js ■ 



= xs«i +(l-a;)s»; 



(C.46) 
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^+ = i + ; R- =R-, 
r = ^ (r + s); f~= r — s. 
Then the Lagrangian L2 gets the form, 



(C.47) 



L2 = EV+V- +x{l-x)Ev+v- + ^£u'^ 



ap 



+ b {f-y + c (^R-y + 2d pf- 



+ 2f pR- + 2ef-R- 



(C.48) 



where the coefficients a, f were defined in (C.41). 
In this case the equation of motion has the form, 



d_ 
dz 



R- = 



f - = 0; 



£ I \ f- + af 
dz, 



0, 



(C.49) 



where 



d _ e ^ 
p+-f- + -R- 
a a 



(C.50) 



Similar to other terms in the Lagrangian which were 
calculated in previous sections, the ^;-dependence of the 
functions R~^{z) and f'^{z) does not affect the corre- 
sponding action. Using the solutions of equations (C.49) 
in the Lagrangian Eq. (C.48) we obtain the action. 



S2 



23 



Z2 














a 




_^ 




a 



T1T2 



' BmiiOi^z') 
{riy+r2~-fj^ + (ff) 





-x{l - 


£ 


{Pl- 


^ 2AZ23 




■-) 




a J 



^ \2 



(Ri)^ + Ri ■ Rr + {RrY 



2i?2 • ^2 + • r\ + i?2 -Ti +Ri ■ 

I 



(C.51) 



where ujq = \/a/£; /^Z2i = Z3 — Z2- 
Eventually, we get. 



where 

A^2(A2;23) = 



W2{z)=N2{Az23)e'^\ 

{E£f 



-i£ujQ 



2-K sin(a;oA2;23) (27rA2;23)^ 



(C.52) 



(C.53) 
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